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We compute the polarized signal from atmospheric molecular oxygen due to Zeeman effect in 
the Earth magnetic field for various sites suitable for CMB measurements such as South Pole, 
Dome C (Antarctica) and Atacama desert (Chile). We present maps of this signal for those 
sites and show their typical elevation and azimuth dependencies. We find a typical circularly 
polarized signal (V Stokes parameter) level of 50 — 300 fiK at 90 GHz when looking at the 
zenith; Atacama site shows the lowest emission while Dome C site presents the lowest gradient 
in polarized brightness temperature (0.3 fj.K/o at 90 GHz). The accuracy and robustness of 
the template are tested with respect to actual knowledge of the Earth magnetic field, its 
variability and atmospheric parameters. 

1 Motivation 

In the last 15 years several experiments managed to successfully observe CMB anisotropies 
from the ground inside the atmospheric windows far from the emission lines of most abundant 
molecules of the atmosphere. For the next generation of ground-based experiments, which will 
be able to map CMB polarization at high resolution and with high sensitivity, an evaluation 
of atmospheric polarized emission is required. As pointed out first by Hanany et aP, the 
main atmospheric contaminant for this kind of measurements is the circularly polarized Zeeman 
emission of molecular oxygen (02)- Such signal can in fact be converted into linear polarization 
through non idealities of the instrument and thus generate spurious Q and U quantities. 
In second instance the circular polarization of CMB itself has a cosmological interest and a 
theoretical effort is undergoing to investigate those mechanisms that can produce such signal 
either at last scattering surface or at later epoch (see for example GiovanniniP Cooray et 
and references in Spinelli et al?). In this perspective any dedicated future measurement from 
the ground aiming to update the current constraint on V Stokes parameter^ has to face the 
presence of atmospheric molecular oxygen as a foreground. 

2 Polarized emission theory 

O2 is the only abundant molecule in the atmosphere having a non-negligible magnetic dipole 
moment since its electronic configuration sets two electrons coupled with parallel spin in the 
highest energy level. This characteristic, together with the presence of the Earth magnetic 
field, allows the Zeeman splitting of its roto-vibrational lines in the millimeter region of the 
electromagnetic spectrum. In particular we notice that O2 has a single intense line around 
v ~ 118.75 GHz and a forest of lines in the range 50 — 70 GHz. Selection rules corresponding 
to transitions Aj = ±1 and Amj = 0, ±1 identify three different type of Zeeman lines (it, o~± 
respectively) having polarization direction and intensity properties varying as a function of the 
angle 9 between the line of sight (los) and the Earth magnetic field. 

The polarized intensity of each line can be evaluated through its coherency matrix p; the total 
coherency matrix Atot is given by the superposition of those for all the different lines weighted 
by two factors taking into account the transition dependence on pressure and temperature 



(C(v, P, T)) and the line broadening and mixing (F(y, u^, Af c ))P 
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Ptrans denotes the transition probability for each line as predicted by quantum mechanics and 
fixes their relative intensities. Taking into account all these relations, the final matrix has the 
following form: 

a 7 f a + b i(a — b)cos6 \ /nN 

+ y— i(a — b) cosv (a + b) cos + csm 9 J 

We notice that the presence of non vanishing, conjugate and purely imaginary off-diagonal terms 
means that radiation is only circularly polarized (V ^ 0,U = 0) with a circular polarization 
stronger when the los is aligned with the Earth magnetic field (V oc cos#). On the other hand, 
diagonal terms are different and thus a small fraction of Q-like linear polarization is produced, 
which is stronger when the los is orthogonal to the Earth magnetic field (Q ex (a + b — c) sin 2 6). 
Nevertheless, the ratio of linear to circular polarization intensity is very low (~ 1CP 4 ) and the 
former can be neg lectedP 

3 Signal computation 

Once the A matrix is computed for a given pressure and temperature configuration at a given 
altitude, we use the tensor radiative transfer approach for isothermal layers to propagate 
the signal through the atmosphere. According to this theory, the signal emerging from one 
isothermal layer is a weighted sum of the physical temperatures of the layer itself. The weight 
function encodes all the information about transfer properties and is a function of the total 
coherency matrix A. 

For the signal computation we used vertical profiles of atmospheric pressure and temperature 
which are publicly available^ and the IGRF-2010 model for the Earth magnetic field. The 
latter describes the geomagnetic field up to an angular scale of ~ 15°, corresponding to a typical 
wavelength of 3000 km along the surface. Smaller scale correction due to magnetized rocks in 
the crust are not accounted for. 

4 Results 

4-1 Frequency templates 

We computed the circularly polarized emission at the zenith for various sites suitable for mm 
astronomy (Chajnantor in Atacama desert, South Pole, Dome C, Testa Grigia) as a function of 
frequency (fig. [TJ. Differences among those sites are mainly due to the difference of altitudes 
between them (lower layers of the atmosphere having an higher pressure are in fact the major 
contributors to the signal) and amplitude and direction of Earth magnetic field, since the signal 
depends on the scalar product between the latter and the los. We notice that both right handed 
(V > 0) and left handed (V < 0) circular polarization are produced and that a sign reversal 
takes place at frequencies v ~ 100, 160 GHz. 

4-2 Angular templates 

We produced maps of the signal at the chosen sites in local alt-azimuthal coordinates at 90 GHz 
(see figure [2] and Spinelli et ai). Both right handed and left handed circular polarization are 




produced varying the los. In some cases a null signal direction is present and corresponds to 
the angular position where magnetic field and los are orthogonal. As expected, the observed 
signal depends on a combination of the atmospheric thickness (elevation scans have a zenith 
secant dependence law) and the magnetic field direction (see fig. [3]). Such dependencies can 
lead the two effects to roughly compensate each other producing a nearly constant signal on a 
large part of the visible sky (see Dome C case in fig. [2]), or to enhance the signal gradient where 
the magnetic field is nearly horizontal, as in Atacama case. 

Typical signal variation at 90 GHz for North-South elevation scans ranges between 70fiK (Dome 
C) and 500/U-fT (Atacama) while for a full 360° azimuthal scan at el = 45° varies between 50fiK 
(Dome C) and 270 fiK (Atacama). We also computed the signal gradient for both type of scans 
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Figure 2: Map of atmospheric O2 circularly polarized signal at v = 90 GHz for Atacama (left) and Dome C 
(right) in Lambert projection centered on the zenith. Black circles denote elevation intervals of 15°; white (red) 
contours denote intervals of 50/j.K in brightness temperature (null signal level). 



in order to give a better description of the signal angular dependence. The maximum value of 



the gradient takes values in a range 0.3 — f>.l[iK/° for elevation scans while for azimuth scans 
it varies between 0.4 — 2.5/xiT/° with minimum and maximum values at Dome-C and Atacama 
respectively. 

Dome C presents the smoother pattern for O2 polarized signal which thus can be minimized for 
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Figure 3: North-South elevation scans (left) and azimuthal scans at constant elevation el = 45° (right) of the 

oxygen polarized signal for different sites at v = 90 GHz. 

any kind of differential scanning strategy. We note however that in case of a realistic application 
the signal computed in our template has to be integrated with the bandwidth efficiency of the 
experiment and convolved with its beam profile. 

5 Accuracy 

The accuracy of the 90GHz templates has been estimated taking into account uncertainties 
and typical variability of the main parameters of the model through Monte Carlo techniques. In 
particular the influence of uncertainties in the magnetic field model direction and magnitude and 
secular variation (SV) have been found negligible for all sites (SV rms < 0.2fiK, < Q.hfxK/y). 
A violent and rapid event like a solar magnetic storm conversely can affect the accuracy of the 
template more significantly, causing a variation of 2 — 9fiK on the level of the V signal. During 
those events in fact the Earth magnetic field strength can vary as much as 1000 nT for no more 
than a few days before quiet conditions are established again. 

We then investigated the uncertainty with respect to the accuracy of temperature and pressure 
profiles using values quoted in literature^ and found it to be below 2fiK for all the sites. Day to 
day variation of temperature and pressure profiles do not affect significantly the signal level while 
seasonal long term variability of atmospheric parameters affects the accuracy of the template 
only for polar sites (3 — 14/x/f), where it concerns the whole air column and not only the lowest 
layers of the atmosphere like, e.g., the Atacama case. 

We note also that the line transition frequencies are known with negligible error bars but the 
Oxygen absorption parameters, which are crucial to compute the line profiles, have error bars 
of at least 5% at 90 GHz. This accuracy applies up to ~ 120 GHz but at higher frequencies, 
where models and values of Oxygen absorption coefficients are not yet validated, results should 
be used with caution. 
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